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Abstract. We present a theoretical profile of the Lyman /3 line of atomic hydrogen perturbed by collisions 
with neutral hydrogen atoms and protons. We use a general unified theory in which the electric dipole 
moment varies during a collision. A collision-induced satellite appears on Lyman j3, correlated to the 
B"B 1 Si — X Eg asymptotically forbidden transition of H2. As a consequence, the appearance of the line 
wing between Lyman a and Lyman (3 is shown to be sensitive to the relative abundance of hydrogen ions 
and neutral atoms, and thereby to provide a temperature diagnostic for stellar atmospheres and laboratory 
plasmas. 

PACS. 32.70.Jz Line shapes, widths, and shifts - 52.25.Qt Emission, absorption, and scattering of ul- 
traviolet radiation - 95.30.Dr Atomic processes and interactions - 95.30.Ky Atomic and molecular data, 
spectra, and spectral parameters - 97.20.Rp Faint blue stars (including blue stragglers), white dwarfs, 
degenerate stars, nuclei of planetary nebulae 



1 Introduction 

In Allard et al. 1999 § we derived a classical path expres- 
sion for a pressure-broadened atomic spectral line shape 
that allows for a radiative electric dipole transition mo- 
ment which depends on the position of the perturbers. 
This factor is not included in the more usual approx- 
imation of Anderson & Talman 1956 Q and Baranger 
1958 We used this theory to study the influence of 

the variation of the dipole moment on the satellites present 
in the far wing profiles of the Lyman series lines of atomic 
hydrogen seen in stars and in laboratory plasmas. 

Satellite features at 1600 A and 1405 A in the Lyman a 
wing associated with free-free quasi-molecular transitions 
of H2 and H \ have been observed in ultraviolet (UV) spec- 
tra of certain stars obtained with the International Ultra- 
violet Explorer (IUE) and the Hubble Space Telescope 

(HST) fHH]]. The stars which show L y man a satellites 
are DA white dwarfs, old Horizontal Branch stars of spec- 
tral type A, and the A Bootis stars. The last two have 
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the distinctive property of poor metal content, that is, 
low abundances of elements other than H and He. Satel- 
lites also have been observed in the laboratory spectra of 
laser-produced hydrogen plasmas Bill 



Satellite features in hydrogen lines are not limited to 
Lyman a, which is the only Lyman-series line accessible 
to IUE as well as HST. Observations with HUT (Hopkins 
Ultraviolet Telescope) and with ORFEUS (Orbiting Re- 
trievable Far and Extreme Ultraviolet Spectrograph) of 
Lyman (3 of DA white dwarfs with T e ff close to 20000 K 
have revealed a line shape very different from the expected 
simple Stark broadening, with line satellites near 1078 and 
1060 A fUjQ. The satellites in the red wing of Lyman /3 
are in the 905 to 1187 A spectral region covered by the 
Far Ultraviolet Spectroscopic Explorer (FUSE) launched 
in June 1999. Furthermore, Lyman (3 profiles are also the 
subject of an ongoing study of the far ultraviolet spectrum 
of dense hydrogen plasmas. The strengths of these satel- 
lite features and indeed the entire shape of wings in the 
Lyman series are very sensitive to the degree of ionization 
in the stellar atmosphere and laboratory plasmas, because 
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that determines the relative importance of broadening by 
ion and neutral collisions. 

In a previous work ||l3f we presented theoretical pro- 
files of Lyman j3 perturbed solely by protons. The calcula- 
tions were based on the accurate theoretical H molecular 
potentials of Madsen and Peek |l4|] to describe the inter- 
action between radiator and perturbers and dipole transi- 
tion moments of Ramaker and Peek flTa] . The line profiles 
were included as a source of opacity in model atmospheres 
for hot white dwarfs, and the predicted spectra compared 
very well with the observed ORFEUS spectra Q . 

The very recent ab initio calculations of Drira 1999 jl6| 
of electronic transition moments for excited states of the 
H2 molecule and very accurate molecular potentials of 
Schmclchcr now allow us to compute Lyman (3 pro- 

files simultaneously perturbed by neutral atomic hydrogen 
and by protons. The aim of our paper is to point out a 
collision- induced satellite correlated to the B"B 1 27+ — 
X asymptotically forbidden transition of H2. We show 
that the shape of the wing in the region between Lyman (3 
and Lyman a is particularly sensitive to the relative abun- 
dance of the neutral and ion perturbers responsible for the 
broadening of the lines. 



the set of electronic coordinates (position and spin) plus 
spin coordinates of the nuclei, while R denotes the set 
of position coordinates of the nuclei. We assume that the 
radiating atom is immersed in a perturber bath in thermal 
equilibrium. The density matrix p is 



,-/9H 



P ; 



where (3 is the inverse temperature (1/kT). 
We use the notation 

(0)=Trp{), 



(5) 



(6) 



2 Theory 

The classical path theory for the shape of pressure-broadened 
atomic spectral lines which takes into account the varia- 
tion of the electric dipole moment during a collision is 
only briefly outlined here. The theory has been described 
in detail in Our approach is based on the quantum 
theory of spectral line shapes of Baranger |||] developed 
in an adiabatic representation to include the degeneracy 
of atomic levels ppaM|. 



where Tr denotes the trace operation. 

The adiabatic or Born-Oppenheimer representation com- 
prises expanding states of the gas in terms of electronic 
states Xe{x', R) corresponding to frozen nuclear configura- 
tions. In the Schrodinger equation 

(T elec + V(x, R)) Xe (x; R) = U elec {R) Xe (x; R) (7) 

= E e {R) Xe (x-R) . (8) 

R appears as a parameter, and the eigenenergies E e (R) 
play the role of potential energies for the nuclei. Any total 
wave function ^(x, R) can be expanded as 



*(x,R) = J2MR)Xe(x;R) 



(9) 



2.1 General expression for the spectrum in an 
adiabatic representation 

The spectrum I(lo) can be written as the Fourier trans- 
form of the dipole autocorrelation function ^(s), 



Here, 



I(u) = -Re [ ${s)e- luJ8 ds 
n Jo 



<P(s) = TrpDte^De - ^ 
= <Dt(0)D( s )> 



(1) 



(2) 



As the nuclei get far from each other, which we denote 
by R — > 00 , the electronic energies E e (R) tend to asymp- 
totic values E^° which are sums of individual atomic ener- 
gies. Since atomic states are usually degenerate, there are 
in general several different energy surfaces which tend to 
a same asymptotic energy as R ~ > 00 . We will consider 
specifically a single radiating atom, the radiator, immersed 
in a gas of optically inactive atoms, the perturbers. For a 
transition a = («,/) from initial state i to final state /, 
we have i?-dependent frequencies 

oVe(-R) = (E e ,(R) - E e (R))/h , e e e l , e' e e/10) 

which tend to the isolated radiator frequency 

u a =Lj fi = (Ef -E?)/h (11) 

as the perturbers get sufficiently far from the radiator: 

uVe(-R) — > w/f as R — > 00 , e 6 c.j , e' G E/ . (12) 

Let us introduce projectors P e which select the e th adia- 



(3) batic component of any &(x,R) according to pfl 



is the autocorrelation function of D(s) , the dipole mo- 
ment of the radiator in the Heisenbcrg representation (we 
use bold notation for operators) [^2| . H is the total Hamil- 
tonian 

H = T„ Md +T eiec + y(x,R), (4) 

where T n ud and T e iec are sums of nuclear and electronic 
kinetic energy operators respectively, and V(x, R) is the 
interaction between particles. Here x denotes collectively 



P e #(x,R)=Tp e (R)Xe{x;R) ■ (13) 

We write the dipole moment as a sum over transitions 

D = ^D Q , (14) 



D, 



e.e' 



P e dp p 



(15) 
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In the Heisenberg representation 
D„(i) 



e,e' 



(16) 
(17) 



mean velocity v, we get from 
ten as 

1 



lLpj that g a (s) can be writ- 



es) 



Ee,e' (Q) Me 



e,e' 



+ 00 



2npdp 



da; dp 



The sum e , is over all pairs (e, e') such that u> e >. e (R) —* u) c 
as i? — > oo. Thus D Q connects all pairs of adiabatic states 
whose electronic energy differences become equal to oj a as 
R — > oo. In the absence of perturbers, D Q would be the 
component of D responsible for the radiative transitions 
of frequency u) a . We note that the projection operators 
account for the weighting factors discussed in Ref. pi] ]. 

Introducing the expansion Eq. (|l4|) for D into the ex- 
pression Eq. (H) for $(s), we obtain 



A Jo dtV e'e[r(t)] 



,[r(0)] 
'r(s)} - d e 



where 



<P a ,f3(s) = TrpB^e"" e 



The line shape is then 



a,/3 



(18) 



(19) 
(20) 



(21) 



The terms I a ^(u)) , a/ft represent interference between 
different spectral lines [Q . When these interference terms 
are neglected, we get 



(22) 



(0)]](27) 

The separation of the radiator and perturber is 
r{t) = [p 2 + (x + vt) 2 J 1 / 2 with p the impact parameter 
of the perturber trajectory and x is the position of the 
perturber along its trajectory at time t — 0. The total line 
strength of the transition is ^ e , ^ \d ee ' | 2 - The potential 
energy for a state e is 

V e [r(t)]=E e [r(t)]-E?; 

the difference potential is 

V e , e [r(t)} = V e ,[r(t)]-V e [r(t)] ; 

and we defined a modulated dipolc 0] 

d eel [r(t)}=d ee ,[r(t)]e-% v ^W , 

where we denoted 

4e'( r ) = (Xe(r)|d|xe'(r)) . 

In the above, we neglected the influence of the potentials 
V e (r) and V e i (r) on the perturber trajectories, which re- 
main straight lines. Although we should drop the Boltz- 
mann factor e~^ v "^ for consistency with our straight tra- 
jectory approximation, by keeping it we improve the result 
in the wings. Note that over regions where V e (r) < 0, the 
factor er^ v '^ accounts for bound states of the radiator- 
perturber pair, but in a classical approximation wherein 
the discrete bound states are replaced by a continuum; 
thus any band structure is smeared out. 



(28) 



(29) 



(30) 



(31) 



and 



where 



= <D£(Q)D a ( S )> 



(23) 



(24) 



The time dependence of <P a (s) is determined by D Q (s), 
the part of the dipole moment which, in the absence of 
perturbers, oscillates at the frequency u> a . Let us now de- 
note 

d a (s) =T) a (s)e^ s (25) 

wherein the free evolution e~ luJaS is factored out. 

For an isolated line, such as Lyman a, we have shown 
(Allard et al. [l)) that the normalized line shape J a (Auj), 
in the uncorrelated perturbers approximation, is given by 



J a (Abj) =FT[e" 9 ° (s) ] 



(26) 



In the classical path approximation, where we assume that 
the perturber follows a rectilinear trajectory at a single 



3 Theoretical analysis 
3.1 Formation of line satellites 

Close collisions between a radiating atom and a perturber 
are responsible for transient quasi-molecules which may 
lead to the appearance of satellite features in the wing of 
an atomic line profile. 

When the difference AV(R) between the upper and 
lower interatomic potentials for a given transition goes 
through an extremum, a relatively wider range of inter- 
atomic distances contribute to the same spectral frequency, 
resulting in an enhancement, or satellite, in the line wing. 
The unified theory 



J,|22J predicts that there will be satel- 
lites centered periodically at frequencies corresponding to 
the extrema of the difference potential between the upper 
and lower states, Aw = kAV cxt , (k = 1, 2, . . .) f||||,|f|. 
Here Au) is the frequency difference between the center 
of the unperturbed spectral line and the satellite feature, 
measured for convenience in the same units as the poten- 
tial energy difference. This series of satellites is due to 
many-body interactions. 
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3.2 Diatomic potentials 

The adiabatic interaction of the neutral hydrogen atom 
with a proton or another hydrogen atom is described by 
potential energies V e (R) for each electronic state of the 
or H 2 molecule (R denotes the internuclear distance 
between the radiator and the perturber). For H-H + col- 
lisions we have used the potentials of H^ calculated by 
Madsen and Peek |l4|]. For H-H collisions we have used 
the potentials of H 2 calculated by Schmelcher jL7| . 

In Fig. [I] we have plotted the H 2 potential differences 
AV(R) for the singlet states which contribute to Lyman (3. 
We have used letters here to label the states. B3 is the well 
known B"B 17+ state. At small internuclear separation, 
R, the state correlates with the B" state, the third 
1 27+ state of the Rydberg series. At larger internuclear 
separation the state has an ionic character until R— 19 A 
where the potential energy curve of the H + + H (Is) 2 
state crosses the H(n = 1) +_H(n = 3) energy levels, be- 
cause of an avoid crossing B looses its ionic character. 
Dabrowski and Herzberg (26| predicted the existence of 
the B, the first calculations were done by Kolos p7| , p8| . 
More recent calculations are reported in M , B9| . 

We use B4 and B5 respectively to label the 4 1 E£ and 
5 states, and V and D to label the 3 1 7Z [i and 2 1 II U 
states. 

Each difference potential exhibits at least one extremum 
which, in principle, leads to a correspo ndin g satellite fea- 
ture in the wing of Lyman [3 (see Sect. 3.1). The present 
approach now allows us also to take into account the 
asymptotically forbidden transitions of quasi- molecular hy- 
drogen which dissociate into (Is, 3s) and (Is, 3d) atoms. 
The satellite amplitude depends on the value of the dipole 
moment through the region of the potential extremum re- 
sponsible of the satellite and on the position of this ex- 
tremum. We have shown |p], p^||3fj| that a large enhance- 
ment in the amplitude of a spectral line satellite occurs 
whenever the dipole moment increases through the region 
of internuclear distance where the satellite is formed. 

The potential differences of the B3-X and B4-X tran- 
sitions exhibit double wells. The maximum at 3.0 A of the 
B3-X potential and minimum of the B4-X potential are 
due to an avoided-crossing. 

The most significant characteristic of Fig. [I] is the exis- 
tence of the deep outer well at 6 A of the B3-X potential. 
The ionic interaction decays slowly making the potential 
energy difference very broad compared to the very steep 
wells of the other transitions. This is very important as 
the position of the extremum and the functional depen- 
dence of the potential difference on internuclear separation 
determine the amplitude and shape of the satellites [ElJ . 



moments for Hj and H 2 , calculated as a function of inter- 
nuclear distance respectively by Ramaker and Peek |l5) 
and by Drira fl^| , were used for the transitions contribut- 
ing to Lyman a and Lyman (3. For Lyman /?, the four com- 
ponents which correspond to Is — 3p atomic transitions are 
dipole allowed Jl(|: the two singlet B4-X and D-X tran- 
sitions, and the two triplet 4 3 S+- 3 S+ and 2 3 II g - 3 S+ 
transitions . 

If electronic states i and / of an isolated radiator are 
not connected by the dipole moment operator, that is 
if Dif(R —> oo) = 0, allowed radiative transitions can- 
not occur between these two states. This happens for the 
other transitions which correspond to Is — 3s and Is — 3d 
atomic transitions. Although these transitions should not 
contribute to the unperturbed line profile, Dif(R) may 
differ from zero when a perturber passes close to the radi- 
ator. In this instance radiative transitions are induced by 
collisions, but not at the unperturbed line frequency. 

It often arises that an extremum in the potential dif- 
ference occurs when the final (or initial) potential en- 
ergy curve exhibits an avoided crossing, the corresponding 
wavefunctions exchange their characteristics and the ra- 
diative dipole transition moment varies dramatically with 
R. This is exactly what happens for the B3-X transition. 
To point out the importance of variation of dipole moment 
on the formation of a collision-induced (CI) satellite, we 
have displayed in Fig. || D(R) together with the corre- 
sponding AV(R) for the B3 — X asymptotically forbidden 
transition. The dipole transition is extremely small for the 
isolated radiating atom ( R — »• oo ) but it goes through a 
maximum at the value of R where the avoid-crossing oc- 
curs and remains quite important at the internuclear dis- 
tance where the potential difference of the outer well goes 
through a minimum. 

In such a case we expect a contribution from this tran- 
sition in the wing and the formation of a CI line satellite, 
if it is not smeared out by larger dipole-allowed contribu- 
tions. 



4 Lyman (3 profile 

The Lyman profiles and satellites shown here are calcu- 
lated at the low densities met in the atmospheres of stars. 
The typical particle densities from 10 15 to 10 17 cm -3 al- 
lows us to use an expansion of the autocorrelation function 
in powers of density as described in [|lj|3l| . Line profiles 
are normalized so that over lo they integrate to 1. 



4.1 Collisional profile perturbed by neutral H 



3.3 Electronic transition dipole moments 

The dipole moment taken between the initial and final 
states of a radiative transition determines the transition 
probability, but for two atoms in collision, the moment 
depends on their separation. This modifies relative contri- 
butions to the profile along the collision trajectory. Dipole 



We will consider the two following mechanisms which con- 
tribute to the Lyman (3 wing. 

— The far wing of allowed dipole lines, due to the free- free 
transition in a pair of colliding atoms. 

— The collision-induced absorption due to the free-free 
transition involving the transient dipole moment ex- 
isting during a binary collision. 



[-H collisions 
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Fig. 1. Difference potential energies expressed in cm 1 for the singlet states of H2 contributing to Lyman beta. 



The line profile calculation shown Fig. [| has been done 
at a temperature of 10000 K for a perturber density of 
10 16 cm~ 3 of neutral hydrogen. The only line feature is a 
broad CI satellite situated at 1150 A in the far wing, due 
to the B3-X dipole forbidden transition. Normally such 
an effect would be overshadowed by the allowed transition 
wing, but in this case there is no large contribution of the 
dipole allowed transitions in this region, as can be easily 
predicted by the examination of Fig. [j]. The extrema of 
the allowed B4-X and D-X transitions occur for very short 
distances, and are much smaller compared to the position 
and dept h of the outer well in the B3-X transition (see 



Sect. 3.1 and 3.2) 



The collision-induced absorption depends on the inter- 
nuclear separation and produces very broad spectral lines 
with a characteristic width of the order of the inverse of 
the duration of the close collision. It is strongly dependent 
on 



— the temperature 

— the amplitude of the dipole for R n 
difference presents a minimum. 



when the potential 



This emphasizes the importance of the accuracy of both 
the potential energies and the dipole moments for the line 
shape calculations. 

Oscillatory structures appear near the 1150 A satellite 
as they appear between the line and the 1600 A satellite 
in both theory and experiment (!]] . These oscillations were 
predicted by Royer p2| and Sando et al. |3^j34|. 



4.2 Simultaneous perturbations by H and H 

The complete Lyman (3 profile perturbed by collisions with 
neutral hydrogen and protons is shown in Fig |^. We notice 
that the collision-induced H-H satellite is much broader 
than the allowed H-H + satellite since the dipole moment 
differs from zero on ly o ver a short range of internuclear 
distances (see Sect. |4.l|) . The CI satellite of Lyman (3 is 
quite far from the unperturbed Lyman [3 line center, actu- 
ally closer to the Lyman a line. It is therefore necessary to 
take into proper account the total contribution of both the 
Lyman a and Lyman [3 wings of H perturbed simultane- 
ously by neutrals and protons and to study the variation 



collisions 




Fig. 2. Difference potential energy AV in cm 1 and the corresponding D(R) in atomic units for the dipole forbidden B3-X 
transition which gives the CI satellite of Lyman f3. 



of this part of the Lyman series with the relative density 
of ionized and neutral atoms. 

In jjj we evaluated both the Lyman a and Lyman (3 
wings of H perturbed by protons. However, we neglected 
interference terms between the two lines. Equation (51) 
of §, which gives the profile for a pair of lines such as 
Lyman a and (3, is 

I(cj)=4>W e n '«<°> J a (u-w a ) 

+ <t>fe n ^)j fj (u- U}0 ). (32) 

The perturbed line strength (j^e n ^ a ^ differs from the 
free line strength (f>a by the factor e n ^° ^ . This density- 
dependent factor expresses the fact that the total intensity 
radiated increases or decreases when the dipole moment 
is increased or decreased, on average, by the proximity of 
perturbers. Because of the low densities we consider, we 
have neglected this factor here. 

We show the sum of the profiles of Lyman a and Ly- 
man in Fig. ||. We can see that a ratio of 5 between 
the neutral and proton density is enough to make the CI 



satellite appear in the far wing. The CI satellite appear- 
ance is then very sensitive to the degree of ionization and 
may be used as a temperature diagnostic. 



5 Conclusions 

In the case of Lyman a and the H-H + Lyman (3 satellites, 
the potential shape played a dominant role in the large 
difference in the broadening of the quasi-molecular fea- 
tures J3(|. The width of the collision- induced absorption 
is determined, for the most part, by the short range over 
which the corresponding transition dipole moment is sig- 
nificant. In the CI satellite, it is the dipole moment which 
is very important and which is responsible of the shape of 
the satellite, the observation of such a satellite would be 
a test of the accuracy of the dipole moment calculation. 
Satellites due to allowed and forbidden transitions depend 
linearly on density. The CI satellite is very sensitive to the 
temperature of the absorption, and it may also be used as 
a diagnostic tool for temperature. We emphasize that the 
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Fig. 3. Line profile of Lyman /3 perturbed by neutral hydrogen and protons. The dashed line ( ) shows only the contribution 

from neutral perturbers. 



effect of finite collision duration does play a role in the 
shape of the far wing. The present calculations are done 
in an adiabatic approximation using a rectilinear trajec- 
tory. This should affect slightly the shape of the satellite, 
although no great error is expected. We are developing 
methods to include trajectory effects in the evaluation of 
the line shape. 
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